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Explicit formulae of the equations in the generalized Galileon models are given. We also develop 

the formulation of the reconstruction. By using the formulation, we can explicitly construct an action 

which reproduces an arbitrary development of the expansion of the universe. We show that we can 

completely and explicitly separate the action to the part relevant for the expansion and the irrelevant 

part. The irrelevant part are related with the stability of the reconstructed solution and we can 

further separate the part to the part relevant for the stability and the part irrelevant for the stability. 

The conditions how the reconstructed solution becomes stable and therefore it becomes an attractor 

solution are also given. Working in the static and spherically symmetric space-time, we investigate 

how the Vainshtein mechanism works in the generalized Galileon model and the correction to the 

Newton law becomes small. It is also shown that any spherically symmetric and static geometry 

can be realized by properly choosing the form of the action, which may tell that the solution could 

C , have fourth hair corresponding to the scalar field. We again separate the action to the part relevant 

3 ■ for the reconstruction for the spherically symmetric and static geometry and the irrelevant part. 

We show that by choosing the relevant and irrelevant parts appropriately, we can obtain an action 

C^ I which admits both the solution corresponding to an arbitrarily given spherically symmetric and 

04 , static geometry and the solution an arbitrarily given expansion history of the universe. 
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INTRODUCTION 



The observation of the type la supernovae at the end of the last century tells that the expansion of the present 
universe is accelerating [ll-Q- In order that the accelerating expansion could occur in the Einstein gravity, we need 
J> ] cosmological fluid with the negative pressure and we call the fluid as dark energy (for review, see j^-Q ) ■ The simplest 
fSJ ■ model of dark energy is the cosmological term in the Einstein gravity, which is called ACDM model (CDM is cold 
\G ] dark matter). The ACDM model, however, suffers the so-called fine-tuning problem and/or coincidence problem. In 
O/ . order to avoid these problems, many kinds of dynamical models have been proposed. Especially there are dynamical 
models using the scalar field(s) like quintessence model |7Hll|. fc-essence models |12h14{ . or ghost condensation models 
CO . I13,ll3l- Such scalar models, however, generate large correction to the Newton law in general by the propagation of the 
^D ' scalar field. In order to make the correction decrease, so-called Chameleon mechanism has been proposed in |19l.[2G|. 
^^ ' In the Chameleon mechanism, the mass of the scalar field becomes large when the mass density is large and therefore 
the force propagated by the scalar field becomes very short range, so that the correction can become very small not 
to confiict with the observational or experimental results. Another mechanism to suppress the contribution from the 
scalar field is the Vainshtein mechanism [2l| , where the scalar field is suppressed by the non- linear structure of the 
scalar field equation. Originally the Vainshtein mechanism was a mechanism for the decoupling of the longitudinal 
C^ , mode in the massive gravity. After that, it was found that similar mechanism works [25|, |26| for the bending mode of 
the DGP model [22h24| |. Then the scalar field models where the Vainshtein mechanism works have been proposed. 
The actions of the original models have a symmetry called Galilean symmetry and hence the scalar field is called as 
the Galileon field |27H29| . The equation of motion for the Galileon field does not include the derivative higher than 
two, which may prevent the existence of the ghost although the condition is not necessary nor sufficient condition for 
no ghost. We also note that the structure of the equation docs not have a direct relation with the Galilean symmetry. 
Until now there are a number of activities related with the Galileon scalar field [30l - [33 | . 

In generalized models of the Galileon field, the field equations of motion is very complicated but it could be useful 
for future researches if the explicit forms of the equations arc given. In this paper, we explain how the field equations 
of motion do not include the derivatives higher than two. As pointed in [31|, the Levi-Civita symbol e^'^P" plays the 
crucial role in the structure. In J32l |. it has been shown that the actions obtained in [3l| are equivalent to those in 
[33|. When we consider the curved space-time, in order to preserve the structure of the field equation, there appear 
the correction terms including the curvatures in the action, which also guarantees that the Einstein equations do 
not include the derivatives higher than two, either. We also develop the formulation of the reconstruction, which 
tells the explicit form of the action reproducing an arbitrary development of the expansion of the universe. We show 
the conditions how the reconstructed solution becomes stable and thus it becomes an attractor solution. We also 
investigate how the Vainshtein mechanism works in the curved space-time for the generalized Galileon model. We 



also show that any spherically symmetric and static geometry can be realized by properly choosing the form of the 
action, which may tell that the solution could have fourth hair corresponding to the scalar field. 

We give explicit forms of the equations as far as we can although they looks very complicated. We believe the 
explicit formulae could be necessary for later applications. 

In the next section, we give a general formulation of the Galileon field and show how the derivatives higher than 
two do not appear in the field equations and the Einstein equation. In Section Hill we consider the FRW dynamics 
and give formulae for the reconstruction. The stability of the reconstructed solution is also investigated. In Section 
IIVI we examine the spherically symmetric and static space-time and investigate how the Vainshtein mechanism works 
for the generalized Galileon model. We also realize any spherically symmetric and static geometry in the framework 
of the Galileon model. The last section is devoted to the summary and discussion. 

II. FORMULATION OF GALILEON SCALAR 

In this section, we give a general formulation of the Galileon models by using the Levi-Civita symbol e^^P'^ in flat 
space-time and curved one. 

We first consider the flat space-time. We now find the contributions from the Galileon scalar tt to the equations of 
motion are expressed as [27H29[ 



^1 — ^ — ^ I ^/il/pcr 
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^3 = -3 {(trn)' - trtf} = --e^fe^-P-d,dpT:dsd,i:, 

£i = -2 {(trn)^ - Strntrn^ + 2trli'^\ = -2€^^^t'"'P'' dpd^Trd^dpTrdsd^Ti , 

£5 = — |(trn)**-6(trn)^trtf + 8(trn)trn3 + 3(trn2)^ -etrn^l 

= -^e^P^'eP-^P'^d^d^ndpd.nd^dpndsd^n. (1) 

Here 11^^ = d^dyir and trll = 11^ . We have also used the following formula: 

-jf'^j^&Pjfi-^jfp + rj^'^rj/^PifPif" ~ r]°'''r]'^''r]^Pf]^P + ?fr]'^''r]''P7j^P - Tj'^'rj'^PTj^pTj^'' + 7j°"'r]'^P'n'-'''Tj^P 

j^^ap Pa jfj, Si> _ ap Pa ju Sfi _|_ ap Pp ^v Sa _ ap pp. jcr Sf ^ ap pu -ya 5p _ ap Pu yp Sa 

-jf^rfPrfi^rfp + r]"''r]^''r]^Pr]^'' - r]°"'r]'^''r]^P'n^P + 7f''r]'^''r]''P'n^P - Tj'"'r]'^Pr]^''r]^'' + tj^-'tj^ Pr]'-'''r]^P . (2) 

Then it is easy to find the corresponding Lagrangian densities as follows 

Cl = TT, 

£3 = ^e^J^eP-'P^d^yTTdpTTdsd^TT ^ 9'^7r9^7rn7r - dPird^Trd^.d^TT , 

£4 = ^ef''^eP''P''dp7rd,,nd^dpnd5daTT 

= i (dpndf'Tr {anf - d^.Trd^'Tid" dpT^dPd^TT - 29^7r9''7ra^9^7rn7r + 29''7ra^7ra''9^7r9^ap7r) , 

£5 = ^e'^^-'^eP'^P'^do^^d^TrOpd^^dydpirdsd^T: 

= i (dpTTdt'TT [U-Kf - 39p7ra^7rn7r9'^9p7r9''9^7r + 2d^,nd''ndPd^TTd''dpTrd''dan - ad^Trd^TrdpO^TT (Dtt)^ 
+3d''Trd''Trdpd^TTd''dpTTd"d^7T + Qd'^Tid^'TTdPd.^TTdpdpTiUTT - 6d''7rd''Trd''d^7rdpdp7Td''d^TT) . (3) 

In the curved space-time, instead of the above Lagrangian densities ([3]), we may consider the following ones: 

£1 = TT, 






= o [^^.^r^^'^T (Dtt)^ - a^Tra^TrVSpTrVS^Tr - 2^^'^l^''^lV ^.dyTiUTi + 2d''TTd''TTVPd^TrV ^^dpiT 






+3d''Trd''TrW pd^nW'dpTTWPd^Tr + 6d''TTd''TrWPd^TrW pdpTraTr - 69''7r9'^7rV'"a^7rV^ap7rV''9^7r) . (4) 

Here V^ is the covariant derivative and the meaning of the suffix "'^*''" in £g and £4 will be clarified soon. In ^, 
^afi^S^f^^pa ig defined by 

^o^p^s^t^^pa ^ g'^f^gP'^g-ypg^'^ - g°'^' g^^" g^"" g^p + g°'^g^pg^''g^'' ~ g'^^g^Pg^^'g^" + g'^'^g'^'^g^^g^" - g^'^g'^^g^'^g^" 
-9"""9'^^g"""9^'' + a'^'^g'^^g^^g^" - g°""g'^''g''^g^'' + g°"'g'^''g^''g^^ - g°-''g^^g^''g^'' + g'"''g^^g'""g^^ 
+g°"'g^"g''^g^'' - g'^''g^'"g'"'g^^ + g^^g^^g^^g^" - g^'^g^^g^^g^" + g°-''g^''g'^''g^^ - g'^'g^^'g^'^g^" 

cxa Bll ^v So , aa Bti -yp Sv a.a Bv -yp 5 a , aa Bu 'Jll Sp aa Bp -yu Su , cxa Bp ^v 6ll /rx 

-g g 9 g +9 g'^g^'g ~g g g^'g^ + g g'^g'^g^-g g'g^'^g +g g'g^g^- (5) 

which gives 

/375gM^pa ^ gt>''g-rPg^'' - gl^"" g'<" g^P + g'^P g''" g^" ~ g^Pg'^'g^" + g'^^g'^'g^" - g'^^g^'^g^" , 






p.^'e'"'" = 2 [g'^Pg''^ - g-'^g'P) , 
V./e'^""' = (^9'^- (6) 

Since Vpgpu = 0, wc find Va [e"i^'>^eP''P'') = 0. By the variation of tt, instead of ([!]), we obtain 

fi = 1, 

£2 = -2yPd„n 



P" 1 



3„ „ „ „ 3 



epf^"""'^ { -^y^dpn\/sd,7T + -R^^sdpTTdxTT 



ff) = e^^'''eP''P^[~2Wpd,7rW^dpTTWsd„7T + ^R\^^d,ndx7r^sdan 

^(0) ^ ^^PyS^p.pa fA^^d^T,^f,d,7TW^dpTTVsd,7r+'^-R\^pd^TTdx7rV^dp7r\/sd,7r 

+ -daTTdpTiR^pp^d\TiR^g^„dTT: - -daTrdpir (V^R^p^^) dxirVsOaTT - -dai^dpirR^pp^'SJ ^d\v:'^ sd^-n: > . (7) 

Here we have used the Bianchi identity 

= VxR"p^^ + "^pR'^p^x + ^t^R^Xp ; (8) 

and the definition of the Riemann curvature: 

[Vp.^AVp = -R\^,Vx. (9) 

Here Vx is a covariant vector. Note that in £5 and 81 , there appear the derivatives of the Riemann curvature, 
which include the third derivatives. By using the Bianchi identity ([5]), we find 



Then if we consider the fohowing Lagrangian densities: 

AA = ^ef''^e''''P''di3Trd,TrR^p^sdx7rd\ , 

AC5 = le''^''^e^'''P''dc.TTd^,TrRp,p^Wsd^TTdxTTd\, (11) 

o 

by the variation of tt, we have 

Af4 - -^e^^-'^e^'-'P^ {2Vpd,TrRap^sdx7rd\ + 8d.TrRap^sVpdxTrd\ 

+2dp7Td^7TR„p^sOTT + 2dpTTd^Tr'V xRaptsd^n} , 

Af5 = le^Pi^eP'^P- h3Wo.d^7rWsd^7rRp,p^dMrd\+^d^7rRp,p^R\^sdr7Tdx7Td\ 

-12df,nRp^l3^\7 sdaTT"^ adxTrd^n + 2daiTdpiTRpyp^^ ^dxir^ sd^ir - 2dandpTTRp^p^R''g^^drnd^n 
-2daTrdpTrRp^l3^'VsdaTTOTr - 2da_'ndp-K^ xRpvii-i^ sda-nd^Tr] . (12) 

Hence, if we define 

A = 4°^ + A£4 , C5 = /:f + A£5 , (13) 

by the variation of tt, wc acquire 

£4 = fi°^+A£4 

= e/^^e'^""" f-2\/f,d,n\7^dpn\/sd„TT + ^R\^^d,ndx7r\/ sd^n - ^dpTrd.nR^^s^pdxn 

--^ pd^irR^p^sdxTTd^TT - d^nR^p^s'^ pdxnd^TT - -di3Trd„TTR„p^sOTT > , 

f , = ei"^ + Af , 
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+ -da'Kdp^TTR^pP^dx'nR'^S^^drTT - -daTTdpTlR^pp^V ydxTlV sdaTT - -V adp-nV sdcrTTRpyll^dxTTd^TT 

3 3 1 

--dandp.TTRpt./^^R'^g^^^drnd^Tr ~ -daTTdpTrRp^p^S/sdaTTDTT^ . (14) 



Note that in the expressions in ([T4| . there does not appear the derivatives higher than two. 
In the curved space-time, wc may also consider the following Lagrangian densities 

Ce = e^^^'eP-'P^dpTrd.TTR^Spa 

= -AfRp^-^gp^R^d^'TTd^n, 

^ 2dpTrd''nDTTR - A^p^^^^'nV''VP^TR^p - 2d''TTd^'TTV,.dpnR 

+8d''TTdfTiV^dP'KRpp - AdPndfTTUTTRpp + AdPTidpTrVd^TrRp^p^ . (15) 

By the variation tt, we obtain 

Eg = -2e^'^'eP-'P^S/f,d,nR^spa, 

Sj = -3e"f'^'eP-'P^Va,dp7TVpd,7TR^Spa + ^e^^'^'e^^'^P'' dpT:dxTrR\^pR^Spa ■ (16) 



Even in the curved space-time, Si's {i — 1,2,- ■ ■ ,7) do not include the derivatives higher than the second ones^. We 
should also note that 



(17) 



Now we find that the general Lagrangian has the following form 

^ 9G^^^(o) _ 1 ^^^^^^ ^^^^ ^ Gei7r)Ce + Gj{7r)Cr + Gg WA ■ 
Here X = 9p7r9^7r and £§ is four times the Gauss-Bonnet invariant: 

£8 = e"^'^'e^'""'R^f,^,R^spa = 4 (i?^ - 4i?^,i?^'' + iJ^.p^i?^'^"") . (19) 

We now show that even in the generalized action (jlSp , there do not appear the higher derivative (higher than second 
order) in the equations given by the variation of n. In fact, we find 



(18) 



£. 



(G) 



dG2 _ 9G2_ dG2^ 



(20) 



AG) 






2-—l-d0ndo,ndjT:dpn{V„d'^nVsd°'7: - V^S'^TrVa^^Tr) 



f)^G 1 3^G 

- n o t d^TrdpTT{2dandaTT\7 sd^T^ ~ 9"7r9a7rV59^7r) + -—-^dandsTrd^ndpTT 
onoX 2 on-' 

dG. 



(y adand-yndpirV sd^n + SdairV ad-yndpTrV sd^n — AdaTrV" dfTrdpTrV sda^ — nTrdyirdpTTy sdc,n) 

- I UaTTCl^nX/ sOpTT - -dyirdpTTS/ sdcrn j - - 
dGs „„ „ „ „ fl „ 3 



dX 
dGs 



daTTdyirV sdp-K - -djirdpTrV sda^ - -GsVjdpTrVsdan 



AG,o) 



dX 



d°'TTdjTrdpTTR^^fdi3TT ~ jGsR^g^^'dandpTT 

da ^^^ d'^irdpirdy-K^J ^fdpT:\J sdaTT - " "" 9^7rV-yOp7rVa-c^cr7r(D7r(9^7r + 2d°'-K\J adpn) 



(21) 
(22) 



9X2 



T4 

dX^ 



.d^G. 



-2dp dan\J ^d°'lTdl3TTdy-K\I sda-n + 2 'VpdaTT'V^d°'7TdpTTd^TT'\/sdaTT 



-d, 



"dirdX 
dCh 
'dX 



d^Trdj^TrdiyTrV sda-TT 



'2\J pd^T:\J ^dpT:\J sda^ + -R^pi^^d^TrdxTrWsdaTT - -dpTTdi^nR^^gW pd\TT 
+^"FFr {daTT'^jd°'TTdi3TTd^nRp„s^dxn + d"'nRp^^^dxndpTTd„n\/sdaTT) 



dX^ 
d^G 






+2— — 2-a"7ri?p„^ dxndijTTdunV sdaTT 



dX 
I dGj 
2~dX 



dndX 



d^ndpndi,TTRp^^ dx^ 



dpTidyTi{V pR\ s)dxn 



(23) 



AG) 



Only the last term in ([23|) includes the third derivative of the metric. This term can be canceled by l^E^ , which is 
given by the variation of — 064 £5 7 



Af, 



(G) 



13^5 pupa 



T9A^T^9^7rc'/37r9,,7ri?^5p^ + -^r— □7ra/37r9^7ri?^5p<^ + -^—d^T:\J xd^T^d^-KR^Spa 



^ In addition to the Lagrangian densities in JTSj, the Brans-Dicke type Lagrangian density 

does not generate the field equation and the Einstein equation inckiding the derivatives higher than the second ones. Here /(tt) is an 
arbitrary function of the scalar field vr. 



IOGa 



2 dX 



IdGj 
8 dn 



dpnduiV pR\^s)dxTr - ^^^^dpiidyiiR-fSpa + ^dpCid^TrR^Spa + ^GaV pdyTiR^spa 



6 

(24) 



The third derivative term in ([25)1 is canceled by the fourth term in ([M]) . 
We also have 



,(G,0) 



al3yS pupa 



1 Q^G 



'G d^G \ 



3 9X2 

1 /„ a^G 



eaTT^x 



daTrdi.nV Rdi.TrV jdpTrV sdaTT 



dX^ 



2 r'^c^TraX 



a^gTT 



dirdX 



y vdii'K ] daTrdnTry -ydpTrS/ sdaT^ 



dG. 



7 



-^ ( --VQa^7rV/3a^7rVT,ap7rV55^7r + -R ^^pdp-KdxTrV ^dpT:V gdai: + -dandpirR pp^dxTrR'^^^^drTr 
1„ „„.„„„. \ a^G. 



5aA, 



--da-rrdp-KR pp^V^dxnVsdaTT j + "olTJ^ nR^^p^drndandp-KV^dp-KVsdaTr 

BG 1 (9G 

+di3-^daTrdpTr\7^dpTTR^„g^dxTr - -^j^dandpTT{\/ ,yR^pij^)dxn\/ sd^TT 



(25) 



Only the last term in (|25p includes the third derivative term and this term is canceled by the seventh term in A£^ 



(G) 



AS, 



(G) 



a.^-y5 Aivpa 



- { dx——d^7i 4- "^^^TT j daTrd^TrW pdiyTrRjSpa + --^dxirV xdaTTdpirV ^d,y7TR^Spa 



4 4 o 



15G5 



laGs 



-daTrdpTiV pd„TrR^Spcr + T;-^9aTrdpn(y„R\0^)dxTT'^5daTT 
ISGs^ „ „ „ ^„ „ 3 



1 9G 

-^-^dxTTdaTTdpTrR^,,p^drTrRjSpa + ^G^dpirR^pJdrTTRjSpa 



(26) 



It is straightforward to see that the equations of motion given by variation of tt from G6(7r)£6- Gt{'k)Cj, and G8(7r)£8 
do not include the derivative higher than two because Gq{'k), G7{tt), and G8(7r) do not contain X. 

We also investigate the equation given by the variation of the metric. The equation corresponds to the Einstein 
equation. Since 



sr;. 



'>9 = gg'^'Sgpu, 
1 






5Rnu 



puXa 



^\^y5gap - Vx^ pSgau - ^a^y5gxp + V ^V pSgxy + SgppR'',^^^ - Sg^pR'' 



pXa 



SRp^ = - [V (VpSg^p + VJgpp) - V^Sg^^ - V^V^ (gP^Sgpx)] 

= I [V^VcJg.p + V.V^p - VHgp, - VpV, {gf^Sg^x) - 2R\P^,Sgxp + K^9p^ + R'Jdpp] 
6R = -dg^.R'^'' + V^V^^, - V^ (gf^^dgp,) , 



we find 
Here 

{n2r 



(27) 



= 77^2 ( nQ'^'R - R""' ) + {^2^ + insr" + (UiT" + {Ur.T" + (n^^T" + {n^T" + {n^y . (28) 



2^2 I 2' 



\y'^G,{.,x)^^-^^dp.d'^... 



(29) 



(^3) 



fj.iy 



-Gsin^X) (a'^Tra'^TrDTT - d'^TTdPTrVdpn - d^ndPTrV^'dpTT) 
1 
2 



+ - {-g^'^nTT + d^Trd^" + ^^'^:^'' - g^^'dWdr) (G3(7r,X)5p7r9''7r) - -Vp (G3(7r,X)a''7rc)^7r9''7r) , 



{Hir 



1 „, r9G4(^^^(0) 1 



(n^r 



-.r 



dX 



c\ 



-G4(7r,X)£e 



a^G4(7r,X)^(o) lc)G4(7r,X) 



9X2 



i-4 



8 ax 



Ce \ d'^TTd^n 



ldG4{n,X) 



{-9^7r5'^7r ((Dtt)^ - VapTrVa^Tr) + 2ap7r9''7rV"a^7rV''a'"7r 



+ 2d''wdPTTy''dpTTaTT + 2d''TTdPTTWdpTTaTT - 2d''TTdPTT\/'' d"" TtW pd^TT - 28'' ndPTrW" dl^TTW pd^TT 

^2dPTrd''TTW''dpTrW''d^Tr} 



+ - {-g^'^nn + d'^nd'' + 8^-713" - g'^'d^Trdr) 



IdGiJTT^X) 
2 dX 



{2dpTTdPTTaTr - 2dPTrd''TrWpdaTT} 



ldG4(n,X) 
2 dX 



{dpirdPn (V^a^Va^TT + V^d'^nd^'n - a^7rV^a''7r) 



^a^Tra^Tra'^TrnTr - 2V^a^7ra^7ra^7ra''7r}] 

-^G4{n,X) i^{R''Pd''Trdp7T + R^Pd^Trdpir) - ^i^a^Tra^TT 

-iVpV' (G4(7r,X)a''7ra^7r) - ^D (G4(7r, X)a^7ra''7r) - ^VpV^ (G4(7r,X).g^'^a''7ra''7r) 



11 

4 



1 

4 
1 



+ iG4(7r,X)i?^'^ap7ra''7r- iv^V (G4(7r, X)ap7ra''7r) + ^D {G4i7r,X)g^"'dpndPTT) 



~^dpTTd„R''''''P , 

oj \ dX ^5 -gG5(7r,X)£7 



a^Gs (7r,X)^(0) iaG5(7r,X) 



ax2 



-£ 



ax 



£7 ^ a^Tra^TT 



1 aGs (7r,X) 



{-a^7ra''7r(n7r)3 + Sa'^Tra'^TrDTrVapTrVa^Tr + 6ap7ra''7rn7rV^a^7rV'"a''7r 



(30) 



(31) 



6 ax 

-2a''7ra''7rV''a^7rV'^ap7rV^a^7r - 6ap7ra''7rV'^a''7rV^a^7rV''a^7r + 3a^7ra''7rV''ap7r(n7r)2 
+Zd''ndPTTVPdpTT{U-Kf ~ aa'^Tra^TrVapTrVarTrV^a^TT - 3a''7ra''7rV'^ap7rV'"a^7rV^a<,7r 
-6a''7ra'"7rVpa^7rV^a^7rV^a''7r - 6a'^7ra''7rV'"a''7rVpa<,7rn7r - ea^Tra^'TrVa^TrVpa^TrDTT 
-6a''7ra'"7rV^ap7rV^a''7rn7r + 6a^7ra''7rV''a''7rVpa^7rV^a^7r + ea^Tra'^TrVa^TTVparTrV^a^Tr 
+6a''7ra'"7rV^ap7rV^a^7rV^a^7r + ea^a'^TrVapTrV^arTrV^a^Tr} 

+ i (-.g^'^DTT + a^Tra^ + a^Tra" - g^'^'d'^^dr) (^^^4^^ (3ap7ra^7r(n^)2 
2 l^ 6 (/A 

-aapTra^Trva^Trva^TT - ea^Tra'^TrVpa^TrDTr + ea^Tra'^Trv^apTrv^a^Tr)} 

+Va il '^'^5(7r,X) ^_35 ^^P^Q^gi'^yMaA^ _ 3ap7ra''7rn7ra^7rV''a^7r + 3a„7ra''7rn7ra^7rV''a^7r 

1^ 6 ax \ '" 

+3ap7ra''7rV'^a^7rV'^a^7ra''7r + 3ap7ra''7rV'"a^7rV''a^7ra^7r - 3ap7ra''7rV'"a'^7rV''a^7ra^7r 

_^aA7ra''7ra''7r(n7r)2 + -a^Tra'^Tra'^TrVa^TrVapTr + 3a''7ra''7rVpa^7rV'^a^7ra''7r 

^-ZdP-Kd" Tx\I pd^'KXI'' d^TxdPix - ■idPixd" txXI pd^-nXI'' dP-ixd^ix - 6a''7ra''7rV'"ap7ra''7rV^aa7r 
+6a''7ra''7ra''7rVpa^7rn7r - 3a''7ra'"7rV''ap7rV^a^7ra''7r - 3a''7ra'"7rV''ap7rV^a^7ra''7r 
+3a''7ra'"7rV''ap7rV^a'^7ra^7r) } 

-iG5(7r, X)/:7 (-2a^7ra''7rn7ri? + \dP-nd''T^\IPd'''nRpa + 4ap7ra''7rV'^a'"7ri?''^ 
8 

+4ap7ra''7rV''a'"7ri?^„ + 2d^TxdP'K\I''dp-KR + 2&''KdP'K\I^dp'KR - ^dP'T^dp'nXI'' d^iiRP" 

-Ad^-KdpTvV^'d^'KRP'' - AdPTTdP-TTV pdaT^R"" - AdPTTd^nypdaTTRP-" - AdPTTdaTrVpdf'TTR'"' 

-4dPTrd^TTWpd''TTR'"' + Ad^TTdPTrDTTR"^ + Ad'^TTdPTraTTR"^ - Sd^TTdPTTWadrRJ'"' - Sd^TrdPirW^drRJ'"' 



+ \ (-s^'^Dtt + a-^TT^^ + ^^'^:^'' - g'^'d^Trdr) \ -lG5{n,X)C7i2dpTTdP7TR ~ AdPird^TrRp^) 
2 L ° 

+Va| - iG5(7r,X)/:7 {-2dpTTdPTrd''TTR^'' - 2dpTrdP7:d''TTR^'' + 2dpTTdPTTd^nR'"' 
-d^TTd''TTd''TTR + Ad^'TTdpTrd'^nRP^ + 2dpTTd''nd''TTRf/' + 2dPnd''nd^'TTRj/' - 2dPnd"TTd^nRJp^')} 
+ (VV^ - g'^'V^) l-^G5{n,X)CTi2dpndPTTDTr - 2dPTrd''TrV pd^n) 

+1 V^ V I - Igs (tt, X)/:7 i-SdpTTdPTrV'df'TT + 8dPT:d''TT\/pd''TT + 8dPTTd''TT\/pd''Tr 
2 [_ 8 

-Sd^TTd^TrnTT)} + iv^V^ I -^G5(7r, X)Cy {-^dpTTdPTrV d^n + SOPird^nV pd^n 
2 L " 

+8dPTTd''TTVpd''Tr - Sd^nd'^TTDTr)} - ^V^ |-iG5(7r,X)/:7 (-49p7ra''7rV'^9''7r 

2 L " 

+49''7ra^7rVp9''7r + 4a''7r9''7rVp9^7r-49'^7ra''7rn7r)}- ig^-'V^V^ (-iG5(7r,X)£7 

2 L 8 

(-4ap7ra''7rV'"9^7r + 8c'''7ra'^7rVpa^7r - Ad'^nd'^TTnTr)} 

+2VpV^( - lG5i7r,X)C7l;{d''TTd''TrV''dPTT + d''7rd''7rV''dP7r 
[_ 8 2 

-29'"7ra''7rV^9''7r - 2c'''7r9''7rV'"a''7r + a^7r9''7rV'"9''7r + a'^Tra^TrVa^Tr) 1 (32) 



{ner = la^'^Gein) 



2- 



-4G6(7r) i-{R''Pd''TTdpTr + R^Pd^-ndpT:) + ^Rd'^nd'^TT 

-4VpV'^ {Ge{TT)dPTTdPTT) + 20 {Gfi{-K)dPTid''n) + 2VpV<, {GQ{'K)gP''dP'Kd''-K) 

-2Ge{TT)RP''dpTTdPTT + 2^^" {Ge{TT)dpTTdPTT) - 20 {GQ{'n)g>"'dpTTdPTi) , (33) 



(Hrr = la^'G^iir) 



+G7(7r)£7 {^23^718" tiOttR + AdpTid''TiVPd''TiRp„ + AdpTTdPTiV^'d'''KR''„ 
+4dpTrdPnV''d''TTRf'^ + 2dt'TrdPTr\7''dpTrR + 2d''ndPTrVdpnR - Ad'^TrdpTrVd^TrRP'' 
-Ad^TTdpTTWd^TTRP" - AdPirdf-KV pd^TiR'"' - AdPird^TrV pd^TiRf" - AdP Trd^nV pd^TiR"" 
-AdPTrd^TTWpd''TTR'"' + Ad''TrdPTraTTR''p + Ad^TrdP-KOTrR^'p - 'ddPiTdP-KV ^drR;""" - id'^TTdPnV ^drRJ^" 

-Sa^Tra^TrV^a^TTi?/^'" - 3a''7ra^7rV^a''7ri?// - 2dpTTdPTTOTTRP'' + 2dPTTd''TTVpd„TTR^"' 



- {-gP^On + a-^Tra^ + a^Tra'^ - gP'^d-'ndr) {G7iTT)C7i2dpTTdP7TR - 4a''7ra'"7ri?p^)} 



-WX< G7{tt)C7 {-2dpTTdPTTd''TTR^P - 2dpTTdPTTdPTTR^'' + 2dpTTdPTTd^TTR'"' 

-d^TTdpTTd^TTR + AdPTTdpird^irRP^ + 2dpT:d"TTd''T:RP p^ + 2dPTrd''Trd''TTRJ' J" - 2dPTTd''T:d^TTR„'' p-) ] 

- [VVP - gP^V^) {G7{TT)C7{2dpTldPTlOTT - 2dPTTd''TTVpdaTT)} 

-iv^V''|G7(7r)£7 (-8ap7ra''7rV'"a^7r + 8a''7ra'"7rVpa^7r + 8a''7ra''7rVpa''7r 
-Sa'^Tra^TrDTr)} + ^V^V^ {G7{-k)C7 {-8dpndPTT\/''d''TT + SOPTrd^TrS/ pd^ir 
f 8a''7ra''7rVpa'"7r - 8a'^7ra'^7rn7r)} - iv2 {G7{tt)C7 {-AdpTTdPTT\/Pd''TT 



+AdPTrdpTrVpd''Tr + AdPird'^TrV pdPn - 4a^7ra''7rn7r)} - -g^'-'V^V^ {G7(7r)£7 






(Hg)''" = 8 (V''V"G8(^)) i? - 8g'^" (V'GsItt)) i? - 16 (VpV^Ggl^)) i?"" - 16 (V^VCsItt)) i?'^'' 

+ 16 (V^GsC^)) R'"' + Wgp" (^pV^Gsin)) i?"" - 16 {VpV^Gs{n)) R'^p'^" . (35) 

Thus, we have given the explicit expressions of the field equations and the equations corresponding to the Einstein 
equation for the generalized Galileon scalar models. 

III. FRW DYNAMICS 

In this section, we give formulae for the cosmological reconstruction by considering the FRW dynamics and we 
investigate the stability of the reconstructed solution. Usually, we start from a theory, which is defined by the action, 
and solve equations of motion to define the background dynamics. The reconstruction is the inverse problem, i.e., 
when an arbitrary development of the expansion of the universe is given, we construct the explicit form of the action 
which reproduces the development. For the reconstructed action, the solution describing the development of the 
expansion history is not always stable. In this section, we show the conditions so that the solution can become stable 
or attractor. 

We now take the FRW universe with the flat spatial part. 



ds^ = -dt^ + a{t)^ J2 {dx')\ (36) 



1=1,2,3 



and we assume n depends on only time variable. 

Let us denote the energy density and pressure of the matters by p and p. Then the first FRW equation is given by 

fc=2 

We also write {Hk)ij as {'Hk)ij ~ 'HkO^Sij. Accordingly, the second FRW equation is given by 

1 ^ 

— (3ij2 + 2i/)+^m. = -|. (38) 

k=2 

If p and p are given by the sums of the contribution of the matters with a constant EoS parameter wi , we find 

p = ^p,a-3(i+-'), p = ^«;,p,a-3(i+-0. (39) 

I I 

Here p/'s are constants. 

We now demonstrate the reconstruction of the expansion history of the universe. That is, for arbitrary development 
of the Hubble rate H{t) or scale factor a(<), we determine the functions G2(7r,X), G3(7r,X), G4{tt,X), G5{tt,X), 
Gq{tt), G7{tt), and G8(7r), which gives the development. Since the redefinition of the scalar field it can be absorbed 
into the redefinition of Gi (i = 2, 3, • • • ,8), we may choose 

n^t. (40) 

Then we find 

dG,{7T,X 



X = -l, G,(7r,X) 



We also write Gi(7r, X) [i ~ 2, 3, 4, 5) as follows. 



, etc. (41) 

t.x=-i 



G,(^,X) = ^Gi")(7r)(l + X)". (42) 



n=0 
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In addition, we suppose the scale factor a{t) can be given by an appropriate function g{t) as a{t) — e^^*\ which gives 
H{t) = g{t). Thus the first FRW equation ([57]) has the fohowing form: 



3 



H- - icf (i) - G^^\i) + {SG^W - lGl>\t)^ m + |-6GW(t) - 6Gf\t) + Jcf (t)} gitf 
+ i^-5Gi'\t) + 2Gi'\t) + ^Gf\t)^g{tf - ISGeimtf - 3QGr{t)g{tf + 48Gs{t)g{tf 
i5^p,e-3(i+-''^W. (43) 



2 

I 

On the other hand, the second FRW equation (|38| has the following form: 



^{sH^ + 2H)+lG^^\t) + lG^i^\t) 

~lGi'\t)g{tf - Gi'\t)git) - G'^l\t)g{t) + \Gf {t)g{tf + ^cf (t)g(f) + \Gf\t)g{t) 

+Gi'\t)gitf + lGi'\t)gitf + Gi'\t)gitym - ^cf WffW' " ^G^^'W^W' - ^cf WffWsW 

-6G6(t)g(t)' - 4G6(i)gW - iGe{t)g{t) + 12G7{t)g{tf + 6G7{t)g{tf + l2G7{t)g{t)g{t) 

^32Gsit)gitf - IGGsimt)' " S2Gsit)mm = -\j2 ^(PiC-3(i+-')9(*) . (44) 

Then if we choose cf (t), G^^)(i), G^ (t), G^^^li), Gf\t), G^^\t), Gf (t), G^ (i), G^W, G'^^\t), GeW, Grit), 
and G8(t) so that the FRW equations (|43| and (|44|) can be satisfied, we have the following solution 

H{t)^g{t), 7r = t. (45) 

We should note that G^'^ {t) (i = 2, 3, • • • ), Gf (t) (i = 2, 3, • • • ), G^*^ (t) (i = 3, 4, • • • ), and G^'^ (t) (i = 3, 4, • • • ) are 
irrelevant for the expansion of the universe. Then we can completely and explicitly separate G„ to the parts relevant 
for the expansion and the irrelevant parts. We should note, however, that some of the irrelevant are related with the 
stability of the reconstructed solution as we see in the following. 

Just for the simphcity, we neglect G4 {7t,X), G5 (7r,X), Ge (tt), G7 (i), and Gs (t) in the rest of this section and 
consider the examples of the reconstruction and the (in)stability. Then, we acquire 



^H' = 2^g"in)il + X) + ^g'-'in)-GfU^) 



K 

+ iX-l) l^-^G^^^ (tt) + 3Gi^) (tt) g' (n) ^G^ (n) g' (n) 

-Y^G^^ (tt) (1 + Xf - 2^G^"' (^) n (1 + Xf'^ ^ + G^^^^if^^ 

n=2 n=2 

~9G3iTT,X)H7t^ , (46) 

-^(3H' + 2h) = lg"(7r)(-l + X)-Ag'2(7r) + Gf (tt) 

+ {X + l) {-\Gf (tt) + 3G^^) (tt) g' (tt) - ^Gf (^) g' (n) 



n=2 

By combining the equation of motion for tt with (j47p , we obtain 



+ ^G^"^ (tt) (1 + Xy - G3 (^, X) TT^ ~ 3G3 (^, X) TT^if . (47) 



A- 



-18 ^ G^"^ (7r)n(l + XY'-^H^tt^ - 6 ^ G^^"^ (^)n(l + X)"-!!/,!^ 

00 00 \ 

6^G^")n(l+X)"-i^4^9^G^"^(l+X)"7TM 
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oo oo oo 

+27 E G^"^ (1 + X)"i/2^2 ^ g ^ ^^(«) (;^ ^ x)"iJ^3 _ 3 ^ Q'{n) ^i^x)"Htt^ 

n—0 n—0 n—0 



2 \ ^g'"{n) ^gT\^) + iGf\^)g'{^) + iG^i\i.)g" {i.) - ^Gf )(^)5'(7r) - ^G^ (7r)g"(7r) 



n=2 

2 



-EGfxi+^r 

Tl = 2 J 

\g\i^)9"{^) - \g"\'K) + f 4ff"'(^) - \gT + 3G(^ V).9'(^) + 3G«(7r).g"(7r) 



-Gf{.)g"{.)\ (1 + X) + G^) + E Gf )(1 + X) 



n=2 

H = -V(7r)(-1+X) + Vw-i.^Gf' W 



y (^ + 1) [-\GT\n) + 3Gi'\n)g'{n) - ^g(")(-)3'(-)) 



.2 °o 



^.2 °o 



y J2 Gt\^ + X)- + y E gT\^){^ + xrn' - -' E G^"H-)"(i + xr-^nH 



n=2 



n=0 



n=0 






n=0 



We can eliminate tt by substituting (|^5)) into (P^ . Here 

C oo „ oo \ 

.^ E Gt\-)n{l + Xr-'^' \n' E Gt\-)i^ + ^)"-' 
n=0 ?i=0 / 

COO oo \ 

-6 E 4"' (-)"(! + ^)""'-' + 9 E G^' (-)(1 + Xrn' 

oo oo 

+42 E G^"H7r)n(l + X)"~^HTr^ - 18 E '^i"^ (^)(1 + X)"HTr 

Tl=0 n=0 

oo „ 

-12j2Gi''\n)n{n-l){l+Xr-'H7:'^ + -g"(7r) - Gf ^(Tr) + 6G«(7r)g'W - 9Gf W.g'(7r) 

n=0 '^ 

oo OO 

+2 E G^")(^)n(l + X)"-i - 4E G^"^(^)"(" - 1)(1 + ^)""' , 



(48) 



(49) 



(50) 



We now define variables x and y hy x = tt and y = g' (tt) /_ff . Then for the reconstructed solution (j45p . we find x = 1 
and y = 1. Since 



dx Try dy q" (tt) o y'^ • d7r yx 

- ' -■ -y x rr^^ff, - 



(51) 



dA^ g' (tt) ' dAT 5'2 (tt) ^ g'2 (tt) ' dN g' (tt) ' 

we can examine the perturbation from the reconstructed solution x — y — 1 and 7r = tasx = l + Sx, y — 1 + 6y, and 



12 



TT = t + Stt by using x and y. Eqs. 



and (1491) as follows 
'dx\ / Sx 



\5y\ ^M\dy 



dN 



i6n I 



, Sn . 



Mil Mi2 Afi3^ 
M = \ M21 M22 M23 
M31 M32 M33, 



(52) 



Here 



+ {9"{t) + 3g'\t) K'Gf\t)) (24Gf (t) - 42GW(t) + ISCf (t)) 

+ (-6GW(t) + 9Gf\t)) S^g"{t) + n' Qcf (t) + 3Gi^\t)g'{t) ^cf (t)5'(t) - ci'^ 

+54Gl'-\t)g'^{t) + 5AG^°\t)g'''{t) - l2G'i\t)g'{t) + l8Gf\t)g'{t) + 24G^^\t)g'{t) 



.(2) 



-6 i^-^g'{t)g"it) -G-r (i)g'W + 3G],'>{t)g'\t) - -G^ (t).9''(i) j + 8G 
+2 (^-lg"'(t) - \Gi"\t) + 3Gi'\t)g'{t) + 3Gi'\t)g"{t) - lGf\t)g'{t) ^cf (t)g"(t) 
M12 = A^'g'-\t) [36Gi'\t)g'^{t) + 6G^^\t)g'{t) + 3 [-6Gi^\t) + 9Gf (i)) g'\t) - 54G,f (t).g'2(t) 
-9G(°)(i)g'(i) + 3Gf\t)g\t) + 6g\t) (^g"it) - Ig^^V) + 3G(^)(t).9'(i) - ^G^"'(t).9'(i) 



K^ 



(53) 



(54) 



Mi3 = A^'g'-\t) -18G^^\t)g'\t)-36G^^\t)g\t)g"{t)~6Gi'\t)g\t)-6Gi'\t)g"{t) 



9"it) + '-g'Ht) - n^Gf + '-g'^t)] (-6G«(t) + 9Gf ) + {-6G^^\t) + 9Gf (i)) g"'{t) 



n(o)(^\„'2i 



(0)i 



ri^^)u\r,"i 



+27G'^'(t)g'\t) + MGt'\t)g'{t)g"{t) + 6G;,"'(t)g'(t) + QGt'\t)g"{t) 
M{t) [ ^9"'{t) - ^-Gf\t)+iG'-i\t)g'{t) + idi\t)g"{t) - ^cf (t)5'(i) - lGf\t)g"{t) 



Ir' 

2 

1 

-( 

2 



-6g"(i) ( ^.g"(i) - ;-G^°\t) + 3Gi^\t)g'{t) -G^ (i)g'(i) 



1 



-2 ( -g""{t) -Gf\t) + 3Gi'\t)g'{t) + 6Gi'\t)g"{t) + 3Gi'\t)g"'{t) - -Gi'\t)g'{t) - 9G]^>{t)g"{t) 



^(0), 



M21 



6 

K 

9"{i) 



6 



+:^5"'W + ^5'(%"'(0-G,^"' 



^(0), 



g,2(i) 5'-'(0 j 5"(0 + «' ( 2^=* ^^^ + ^G\'\t)g'{t) - -G]^>{t)g'it) - Gi'>{t) 



-g'-\t)(^-n'Gi'\t) + L'G^^\t)yiii, 
M22 = -g'-\t)g"{t) - 3 - g'-Ht) [-K'Gi'\t) + ^Ac^Gf (i)) M12 , 
M23 = -g'it)-' [-^'Gi'\t) + ^K^Gf (t)) 7\/i3 + 3g'~\t)g"{t) - g"'^{t)g"'{t) , 



(55) 

(56) 

(57) 
(58) 



where 



A, = -^ (6GW(i) - 9Gf (t))' + i2Gi'\t)g'{t) - 27Gf\t)g'{t) - 2iGf\t)g'it) + ^g"{t) - cf {t) 



K^ 



and 



-6Gi'\t)g'{t)-8G^^\t), 

M31 = 5'"' (t) , M32 = 5'"' (0 , ^^33 = -g'"' (0 3" (0 • 



(59) 
(60) 
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The eigenvalue equation for the matrix M in ((52|) has the following form: 

A^ + aA^ + ^A + 7 = , 
where 



(61) 



a = -A/n - A/22 + 9'-^ (t) g" (t) , (62) 

/3 = -(A/n + Ah2)9'-^{t)g"it) - (A/23 + Ah3)g'''{t) + A/nA/22 - A/12A/21 , (63) 

7 = (A/iiA/22 - A/i2A/2i)5'"'(t)g"(t) + (A/iiA/23 - A/13A/21 - A/12A//23 + A/i3A/22)5'"' (0 . (64) 

When all the eigenvalues A are negative, the reconstructed solution becomes stable. The Hurwitz theorem tells the 
condition that all the eigenvalues are negative is given by 



(i) a > , (ii) a/3 > 7 , (iii) a/3-/ > 7^ . 
For the simplicity, we further put G3 (t) ~ G^ (t) ~ 0. Thus we find 

Mn = A^'g'-\t) [72G^^\t)g''{t) + 24G^^\t)g'{t) + 24Gi'^ {g'\t) + 3g'\t)) + 24Gi'\t)g\t) 
-^g'{t)9"{t)+8Gi'^ + ^g"'{t)\, 

Mi2^6A^'^g"it), 

Ah^ = -2A^'g'-\t)^g""{t), 

hi 

A/21 = , A/22 = -g'-' (t) g" (i) - 3 , A/23 = ^g'-\i)g"{i) - g'-\t)g"'{t) , 
A/31 = M^2 = g'-^ (t) , A/33 = -g'-^ (i) ,9" it) . 



Here Aq has the following form: 



A, = -24Gf (t)5'(0 + ^9" it) - 8G^'^ 



The eigenvalue equation becomes 

(A/22 - A) {A^ - (A/11 + A/33) A + A/iiA/33 - A/13A/31} = . 
Therefore we obtain 



A = A/o 



1 



A/u + A/33 ± V (A/u + A/33) - 4 (A/iiA/33 - A/13A/31) 



Consequently, the condition of the stability is given by 

(i) A/22 < , (ii) A/11 + A/33 < , (iii) A/11M33 - A/13A/31 > . 
From the first condition (i) in (j70p . we find 

0>-5'-2(i)g"(t)-3, 
and from the second condition (ii), 

> A„ 1 [-g'-'{t)g"{t)Ao + UAG'^^I) g' (t) + 2AG^^\t) + 2AGi^\t)g'-\t)g"it) 
6 „,., ,. ■^■?^ ,., ,1 ,., 2 



+24G ^^(t) - -g"{t) + 8G','>{t)g'-\t) + -g'" {t)g'-\t) 

The third condition (iii) in ((70)) presents 

< A-,'g'-\t) [~g'-\t)g"{t)[l2Gf\t)g'\t) + 2AG^;^\t)g'{t) + 2AGf\g"{t) + ig'\t)) 



(65) 



(66) 

(67) 

(68) 

(69) 

(70) 
(71) 



(72) 
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+2AG^^\t)g\t) - ^9\t)g"{t) + 8G^^\t) + ^9"'{t)\ + 2g'-\t)^r'{t)\ . (73) 

Before going to non-trivial case, in order to check the above formulation could work, we consider the de Sitter 
space-time, where g' (t) is constant. Then we find 

^j^( ^" _°3 l\ j^,^ ^ U4G^^\t)g''it) + 24G^i\t)g'{t) + 24G^^\t)g' jt) + SCf 

V 9'-' it) 9'-' (t) ; ' " -24Gf (i)g'2(i) + ^g'{t)g"{t) - 8G'^\t)g'{t) 

Then by solving the eigenvalue equation (|68p . we find the eigenvalues: 

;^^_3 lUG^^\t)g'Ht) + 24GPit)g'it) + 24G^^\t)g'it) + SG'J^ 
-24G^^\t)g'^{t) + ^g'{t)g"{t)-8G^^\t)g'{t) 

Note that the eigenvector corresponding to the eigenvalue is given by 

(76) 




which tells that the eigenvector with the eigenvalue corresponds to the shift of n or the origin of time and therefore 
the eigenvalue does not corresponds to any instability. Especially if we consider the case that Gg (i) = 0, the 
condition that the last eigenvalue in ([75]) is given by 

- '^'"Wf'"'" -.<0. (77) 

G'i> (t) 

can be satisfied if 

G^^^{t)=Dcxp{-f{t)), (78) 

where D is a constant and /(t) is an arbitrary function satisfying the condition < f'{t) < 3g' (<). We should also 

note that Eq. (|78p can be satisfied if Gj (t) is a constant. 

As a little bit non-trivial example, we study the case that g' (t) is a solution of the ACDM model: 

a (t) = ^ sinhi (6i) , (79) 

where A and b are positive constants. Eq. (|79p gives 

9 26^ 

-booth (bt) , g" (t) = ^ 

3 V ; , y V y 3sinh2(6i) 

Then the first condition (i) in (1701) gives 



g' (t) - ^fe coth {bt) , g" (i) = - ^ . ^„ ,,^ . (80) 



2 cosh^ [bt) 
which is trivially satisfied, and the second condition (ii) presents 



^>7. TTTTT^^, (81) 



" \ 3 cosh(6i) sinh(fet) 3 \^ cosh'(&i)/ 



cosh (bt) , 
tanh(6t)G^'^ + -^ ( —^-- ,, " .,, , ) } ■ (82) 



Here 



■— tanh(6t)G^^^ 4- -^ ^ ^ ^ 

b "■ ' ^ k2 \^sinh2(^i) cosh2(6i)sinh^(fe0 



Ao = -16Gf 6coth(6t) - 1^ . i^^ . - 8G^') . (83) 

oK^ smh (ot) 
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By using the third condition (iii) in ((70|) . wc find 

< A^'^ \64b^coth'^{bt)G'^^ + 16bcoth{bt)G'^^ 



-^ — G^'' + 166coth(6t)G^'^ 
sinh2(6i) ^ ^ ' ^ 



-8G. 



(2) 



16 

^ ( — coth(5t) 



8 cosh(6t) 
3sinh^(6i) 



(84) 



Although the above expressions (|5^ and (|55)) are very comphcated, these conditions can be satisfied. The simplest 
example is given by 



Here G is a constant. By substituting ([55]) into ([51]), ([53]), and ([51]), we find 

4&2 1 



An 



3k^ sinh^ (bt) 
> An '1(^24 - 



8G, 



12 



< AnM 166 G 



cosh^(6t) 
62 



c + ^ 



) coth(6t) - 



sinh'^ (bt) cosh^ (bt) sinh^ (6t) 
863cosh(6i) \ 



(85) 

(86) 
(87) 



3^2y-"""v-/ 3K2sinh3(6i)J ■ 
If G > 0, we find Aq < in (|86| and the inequality ([87]) is satisfied. The inequality ([88]) is also satisfied if 

Then stability can be realized. Note that even for the de Sitter space-time, if G2 (t) and G^ (t) are given by 



(89) 



the stability is realized as discussed after ([78| . 

As a result, in this section, by exploring the FRW dynamics, we have given the explicit formulae for the reconstruc- 
tion and we have investigated the condition of the stability for the reconstructed solution. 



IV. VAINSHTEIN MECHANISM 



In this section, in order to investigate if the Vainshtein mechanism could work, we consider the behavior of the 
Galileon scalar field it in the spherically symmetric space-time, especially in the Schwarzschild background. 

The Schwarzschild space-time has a spherically symmetric and static metric. The general spherically symmetric 
and static space-time has the metric of the following form: 



ds^ = -C2*di2 ^ g2A^^2 



'■{de^ + sui^ 



'■), 



(90) 



where $ and A only depend on the radial coordinate r. In the Schwarzschild background, the equation derived by 
the variation of the Galileon scalar field tt is given by 



1 Sy^G2{7r,X) ^ 1 <57^G3(7r,A)£3 , 1 ^^^^^#^4°' 1 1 



S 



-9 
1 5y 



Sit ^^f- 

— 0G5(7r,X) ^(0) 
■9 dX '-5 



-g Sn -y/— g 5tt 8 \/—g ^ — gG ^{t: ^ X) Cq 

1 <5y^G5(7r,A)£7 , 1 5./^Gq(t:)Cq , 1 6,/^G7{tt)C7 



Sn 



/—g Sn 

1 S^^Gs{n)Cs 



-g 



Sn 



-g 



Sn 



-g 



Sn 



Sn 



= 0, 



(91) 



SA dA dA 

where we have used the notation —— = — dn-^rv-r — r. The explicit forms of fj, and also ('Hk)ui' ^-i'^ given in 

Sn on ayo^n) 

Appendix [B] 

We also note that the Einstein tensor is given by 



Gt 



1 2* '^ r /I 

J.2 fj^j. 



2A 



)], 
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Grr - 


- -4=" 


r dr 






Gee = 


-- r'e-" 


dr^ ' \dr J ' r dr 


d^dK 
dr dr 


IdA 



G^^ = sin^^Gee. (92) 

The Einstein equation ([^5)) tells that in order that the modification of the Schwarzschild geometry due to the Galileon 
scalar could be small, we find 

8 

E(^Om^^O. (93) 

j=2 

We now investigate how the Vainshtein mechanism works for the generalized Galileon scalar model. In the limit of 
r ^ GAf , the Schwarzschild metric 

ds^ ^-(l- ^) dr^ + -^4gm + '■'(^^' + '"^' ^'^'^') ' (^4) 

for general spherically symmetric and static space-time in ((90)) behaves as 

$(r) ~ , A(r) ~ . (95) 

r r 

If the Vainshtein mechanism works, the gravity with the Galileon scalar behaves as the usual Einstein gravity in the 
short distance compared with the cosmological scale and $(r) and A(r) in ([90)) behaves as those in (|95)). 

Before examining the Vainshtein mechanism in the covariant and generalized Galileon model, we review how the 
Vainshtein mechanism works when we only include £4 in (j3]). In this case the energy momentum tensors of the 
Galileon scalar behave as 

(ttO^tt" f7r')4 (7r')4 

T- ^8c4^^+C4^ + ---, T[, ^ ~5c,^—f + ■ ■ ■ . (96) 

If wc write C4 = Ci-r^, we find £4 ^ 0(1)- 

Wc now study the generalized Galileon model in ()18p but just for the simplicity, we only include £2 and £3: 

C^^-^C2X + G3i7:,X)C3. (97) 

We may assume the first term in (|97)) dominates in the short distance but the second term does in the long distance. 
Since the field equation has a form like 



when r <g^ ry, we find 



where 



n"{r)^f{7T{r),7T'{r))^^, (98) 



n'{r) = n'iry) exp ( / ZMllli^V)) ^, A ^ r:'{rv)^ , (99) 



Fir)^l f(-i^')y(^'))dr. (100) 



If we provide 7r'(ry) ~ $^(ry), 7r'(r) becomes small when we choose /(7r(r), 7r'(r)) so that we can have F{r) ^ F{rv)- 
If 7r'(r) ^ 7r'(ry) as required so that the Vainshtein mechanism could work, we may find 7r(r) <C 7r(rv'). Then if 
we can choose Gs^n^X) so that /(7r(r),7r'(r)) > when 7r'(r) < n'^ry) and 7r(r) < 7r(rv'), we may surely obtain 
7r'(r) <^ 7r'(ry), consistently. We now have 

... ,_, _2r2-2^(^')'^'-3G3Vr-2^(7r')V-2|%(7r')V , ^ 

fiTTir),TT'{r)) = 2^^^ '-jr^ ^-^^ ^,^ d7^d2£^J 101 

r2 + 6G37r'r + 14^(7r')3r + 4f^(7r')5r 
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It is easy to find that there exists G3{tt,X) which satisfies the condition that /(7r(r),7r'(r)) > when 7r'(r) < 7r'(rv) 
and 7r(r) < Tr{ry) but G'3(7r,X) cannot be uniquely determined. 

Finally we show that any spherically symmetric and static geometry given by arbitrary $(r) and A(r) can be 
realized by properly choosing G^'s. Since the redefinition of the scalar field can be absorbed into the redefinition of 
Gi's, here we may identify the Galileon scalar field with the radial coordinate, tt = r. Then we find 

9G4(r,e-2A)dAe2*-6A d f dG4{r,c-^^)\ e^"^-^^ 1 aG4(r,e-2A) e^^-^A 

—5- ' 



dX dr r dr \ dX J r 2 dX r^ 

^-G4r,e 2Aw + L^ L + -G4r,e ^A ^+ G4r,e ^^ 

2 dr r 2 dr r i r^ 4 

1 d /9G5(r,e-2A)\ e2*-8^ 7 SGsCr.e-^A) ^A e^*-^^ 3^, ^a^x^Ae^^-^A 



2drV dX J r^ 2 dX dr r^ +4<^5(?',e , ^^ 

15^. _2A,dAe2*-6A ldG5(r,e-2A)c2*-4A ^ 3rfG5(r,e- 



4 dr r^ 4 dr r^ 4 dr r^ 

+ 12G6(r)- 4 7 2G6(r) ^ 2G6(r) ^ 

dr r dr r r^ r^ 

^^ , ^dA e^^-^A ^dA e^^-sA dG7(r) e^^-^A dGrW e^^-^A 

-6G7(r)- 2 — + SOGtW- 2 — + 2^ 2 ^^ :^— 

dr r^ dr r^ dr r"^ dr r^ 

_192^^^£!!l!^^ + ^gdGsW 2!!:^^ _ ^grf^GsW 2!!:^ _ ^g dG8(r)e^^-4^ dA 

^j, j,2 ij^ ^j, j,2 ijj, |jj,2 j,2 ^j, j,2 ^j, 

d^Ggfr) e2*-4A 
+16 ^;r ,2 > (102) 

-L|-le-(l-e-) + ^f 
2k^ [^ r^ r dr 



re + e 



dX \ r^ \dr J dr 

aG4(r,e~^^) d$e-4A 5 aG4(r,e-^^) e'^A c>^G4(r,e-^^) d$ Q-e^ a^G4(r,e"^^) e'^^ 

dx 'di^~ 2 ax ^^ dx^ d7^^ ax2 ^ 

_3g4(^ e-2^) — ^ - -G4(r e-^^)^ + i (^iillf:!^ _ 5G4(r^e^d$ c^ 

2 ' dr r 4 ' ^,2 4 ^,2 g^ ^j,, ^ 
1 9G4(r,e-^^)e-2^ 1 aG4(r,e-^^) e"^^ 7 aG5(r,e~^^) e'^A ^^ d^G5{r,e-^^) e'^^ d<^ 

"2 ax ~^2 dx ^~ 2 ax ^d7 ax2 ^d^ 

3 -2An^£1^_1^^ / .-2A-,^2l!^ 1 aG5(r, e-^^) d$ e^^A 3 qq^(^^^ g-2A-) ^^ g-6A 

-M^l^-48^'l?e--. ' ' ' ' (103) 



As in (|42]) , if we wrote 



00 
G2(^,X) = ^G^")(7r) (l-e2A('-='^)x)" . (104) 



n=0 
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wc find 

G2(r,e-2A)^Gf(r). 
Then Eqs. (|102p and (|103p can be rewritten as 

Gf (r) = F(")(r) 



9X -^2 ^"^^ 



(105) 



2e" 



-2* 



ZK'^r'' ar 



-2A 



)] 



^2*-2A 



(G3(r,e- 



,aG4(r,e-^A)dAe' 



-,2$-6A 



ax 



dr 



2 dr 
_d_ /aG4(r^-2A)N^ e2*-6A 



-2A)e-2A) 



19G4(r-,e-2A)e2*-6A 
2' 



G.(.,e-)^2!!:^ + 1 dG,(^^e-- ^ 1^^^^^^_,, 



ax 

„27r-4A 



1 



dr 



v 



1 d / 9G5(r,e-^^) 

'2d7 1 ax 



ic J A p2*-6A 



r 2 dr r 4 

2 



G4(r,e 



-2A\ 



-4$-2A 



dX 



dr 



1 dG5(r, e-^A) j,2*-4A 3 ^^^(^^ g-2A) g: 



o JA .2*-4A 

-2h\ „2*-6A 



-12G6(r 
-6G7(r) 



dA e^^-4A 

dr r 
dAe^*2!^ 
dr r^ 



dr 
dG6(r) e^^-^A 

dr r 

dA e2*-6A 



r^ 4 

„2$-4A 



dr 



r^ 

-4*-2A 



2G6(r)- 



-2G6(r)- 



30G7(r 



dr 



„dG7(r)e2*-4A ^G7(r)e2*-6A 

2 ; 6- 



dr 



dr 



dGgfr) e2*-2A rf$ dG8(r)e2*-2AdA 
-192 — P-^ + 16 — ^ ^ - 16 



d'GsWe 



2<I>-2A 



dGsir)e 



2*-4A 



dA 



-16 



dr 
d^Gsir) e 



dr 



dr 



dr 



dr^ 



dr 



dr 



2*-4A 



dr^ 



G^^Hr) = F^^Hr) 



— ^„2A 



F(o)(r) 



1 
2^ 



-:^e-(i 



,-2A 



+ ^e2^G3(r,e-2A) 

dG3{r,e-^^ ) fe 
dX 

-G3(r,e-2^) 



-3A 



) + 



2d$ 
r dr 



dr 



fr^e*-^^ 



-4A 



dA\ 
d7/ 



-3A 



r^ \ dr 



fr^e*-^ 



))+«-"^ 



r^ 



\dr ^ 



r^e*-^^ 



2e 



-4A 



-TT 



dr 



fr^e*-^^ 



dA 
dr 

-2A\„--2A 



l)(G3(.e-K-)-i»--A(G3(.-.e-)) 



.5G4(r,e-2A)d$e 



-4A 



5aG4(r,e- 



,a2G4(r,e-2A)rf$e 



-6A 



ax dr r 2 

a2G4(r,e-2A)e-6A 3 



ax r2 

2A.rf^e-2A 



ax2 



G4(r,e-^A) 



dr 



-G4(r,e 



aX^ dr r 

-2A)C-^-^ ^ l (G4(r,e-2-^) 



-4A 



aG4(r,e-2A) ^^g 

ax dr r 

7aG5(r,e-2A)c-6Arf$ 

2 



1 aG4(r, e-2A) e-2A 1 ^^^(r, g-^A) g-^A 
2' 



ax 



ax 



ax 



dr 



a2G5(r,e-2A)g-8A^^ 3 _2A,d$e-2A 



15^ , -2ANd*e-'^'^ 



aX2 r2 

iaG5(r,e-2A)d$e 



dr 

-4A 



-12G6(r 



d<i> e 



-2A 



dr 



6G6(r) 



ax dr r 

e-2A 2G6(r) 



, JGs{r) 1 d$ ^^ dG8(r) 1 d* 
-16 ; TT— 48 ; T-;— e 



r^ 

-2A 



6G7(r) 



dr r^ 
3 aG5(r,e-2A) d^e-^^ 

2 ax 

d$ e^2A 



dr 



30G7(r) 



gj^ j,2 

d$ e 



-4A 



dr 



dr r^ dr 



dr r^ dr 



(106) 



(107) 
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Then arbitrary geometry given by $(r) and A(r) can be realized for arbitrary functions Gi (i = 3, 4, • • • , 8) by choosing 

6*2 (r) and (5*2 (r) as in (|106p and (|107p . This may also tell that the spherical symmetric solution could have fourth 
hair corresponding to the scalar field in addition to the usual three hairs corresponding to mass, angular momentum, 
and electric charge. Then only Gj and Gj in G2(7r,X) are relevant for the reconstruction of the spherically 
symmetric and static solution and G2 (tt), n = 2, 3, 4, • • • are irrelevant for the reconstruction although they may be 
related with the stability of the reconstructed solution. 

As we have succeeded in the reconstruction of arbitrary spherically symmetric and static geometry, we may show 

(n) 

the reconstruction is compatible with the cosmological reconstruction in the last section. Let us assume G, , n = 
0, 1, • • • ,Ni in (|42t are relevant for the cosmological reconstruction and G^" , n = Ni + 1, A'^^; + 2, • • • are irrelevant. 
We also assume G,'"' , n = 0, 1, • • • iV^ as in pM| 

G,(7r, X) = Y1 ^."^ (^) (1 - c^^^'^^Xy . (108) 

n=0 

~ (n) 

are relevant for the reconstruction of spherically symmetric and static geometry and G] , n = Ni + 1, Ni + 2, ■ ■ ■ are 
irrelevant. By expanding (g^ as a power series of 1 - e^^^-^^X by using 1 - AT = 1 + e-^Al-A) _ e-2A(0) (^i _ e2A(0)x) , 
we find 

2~k 



G,(7r,X) = ^^Gi")W„Gfc(l+e-2AW)) (-l)^e-2'=A(*) (l - e^^^^XJ 

00 00 , _ 

^(-l)"e-2»A(^)^Gf)WfcG„(l+e-2^(«) '"" (l - e^^'^x)" . (109) 



rGr^(^)fcG„(l+e-^^^^);' " 

n=0 k 

By comparing the expression in (|109p and (|108p . we find 



ex: 7 _ 

G(")(7r) = (_l)"e-2"AW Y. fc^" (1 +e-2A(*)j ' "g^ (tt) . (110) 



Then when n < Ni, we find 



G'f\^) - (-l)"e-2"A(0) ^ ,.c„ (^1 +e-2A(0)j ■ "Gf^TT) 

°^ 1 _ 

= (-l)"e-^"^(^) ^ ,G„(l+e-^^(«)'""Gp)(^), (111) 



k=N i + l 



and when ri > Ni 



G|"^W = (-l)"e-2"''^*) 5] feG„(l + e-2A(^)j Gf)(^). (112) 

k=Ni + l 



k-r 



:(")/' 



For a set given by G^ (tt), n = 0, 1, • • • A^i and G,^ (7'')i f^ = 0, 1, • • • AT;, we can always choose G^ (tt), n = Ni, N + 
1, • • • , A^, + TVj to satisfy piU)) (we may set G*"^(7r) ^ 0, n ^ N + N,N, + N + !,■■■ and G|"'^(7r) = 0, n == 
A^^i + 1, A"i + 2, • • • ). Therefore we can obtain an action generating both an arbitrary spherically symmetric and static 
geometry and an arbitrary expansion history of the universe, simultaneously. 

V. SUMMARY 

In this paper, we gave explicit formulae of the equations in the generalized Galileon models. Even in generalized 
Galileon models, the derivatives higher than two do not appear in the field equation of motion nor the Einstein 
equation. These structures are clearly given by using the Levi-Civita symbol t^'^P'^ . We have also developed the 
formulation of the reconstruction. As a consequence, we can construct an explicit action, which reproduces an 
arbitrary development of the expansion of the universe. We should note that the functions G„'s can be completely 
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and explicitly separate to the parts relevant for the expansion and the irrelevant parts. Some of the irrelevant parts are 
related with the stability of the reconstructed solution and we can also identify which parts in C?„'s are relevant for the 
stability. Then we have succeeded to show the conditions how the reconstructed solution becomes stable and therefore 
it becomes an attractor solution. Working in the spherically symmetric and static space-time, we have investigated 
how the Vainshtein mechanism works. It has been also shown that arbitrary spherically symmetric and static geometry 
can be realized by properly choosing Gi's, which may tell that the solution could have fourth hair corresponding to 
the scalar field. We again identified the parts in in G2{tt,X) relevant for the reconstruction. We should note that by 
choosing Gi (tt, X) appropriately, we can obtain an action which has both the solution corresponding to an arbitrarily 
given spherically symmetric and static geometry and the solution an arbitrarily given expansion history of the universe. 
The generalized Galileon model contains many functions denoted by Gi, which cannot be determined by the ex- 
pansion history of the universe, the stability of the reconstructed solution, nor the condition that the Vainshtein 
mechanism works. These function could be restricted if we consider the developments of the several kinds of fluctua- 
tions by density perturbation, etc. The analysis by the perturbation could be very complicated in the models but it 
might become possible by the development of the technologies. 
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Appendix A: Explicit forms of S^ and {'Hk)i.iv in the FRW universe 

In this appendix, by assuming the FRW universe ([5^ . we give exphcit forms of f^ in ([7]), p^ . and p^ . and (Hfc)^i/ 
in (Uni). 

In the FRW universe ([5S|) . we have 



r* 



a^HSi^ , F 



R 



itjt 



]t 



r: 



H5\, , 



H + H^] a^Srj , R^Jkl = a^H^ {5,k5ij - 5u5k,) 



(Al) 



Then we find 



fi = 1, 

£3 = -ISHnn - 27ij2^2 - 9H'k'^ , 

4°^ = 36H^T:^n + SeH^n^ + SOHHn^ + (third derivative terms) , 

A£i = mH^TT^TT + ISH^-k^ + GHH-k^ + (third derivative terms) , 
£. 



(0) _ 



Afs 



-20H-^7r'^n - 15i/**Tr* - 2lH''Hii^ + (third derivative terms) , 
-SOH^n^n H'^tt^ H^Hn'^ + (third derivative terms) , 



£q = 24H^n + 72H'^TT + 48HHTT, 

£7 = -72H^irn - 108HV - lOSH^Hir^ , 



(A2) 



We also obtain 

('^2)00 

i'H2)ij 

('^3)00 



1 dG,{7r,X) 

-G2{TT,X)a^5ij , 



(A3) 
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("^4)00 

('^5)00 



^^ ^ --G3(7r,XV^ 



a 5i 



-b ^— H TT -6 -— ^ H IT +-G4(7r,AjiJ tt , 

■ 3 9G4(^,X) , , _ d_ ( dG,in,X) \ , dG,{n,X) ^ , ^ dG,{n,X) , 
2 dX dt \ dX J dX dX 

f ^G4(^, XjH^n^ + i^^iiZ^ifTT^ + 1g4(^, XjHn^ + Gi{n, XjHnn 

JG.JTT^X) d^G,{7T,X) 15 3 3 

-5 ^^ H n + ^^2 H ^ +—G5{n,X)H n , 



a Si. 



^3,5 l^pC.j.X) 2.^ 



2 dt V 9^ 



9G5(7r,X) ^^.5 , 5 9G5(7r,X) ^^,.,.. 

dx ^^" + 2 ax ^ " " 



-^G5(^,X)iJ^^^- ^ ^^^l^^'^^ iJ^TT^- ^G5(7r,X)gg^^- ^G5(7r,X)g^A 



a^% , 



(^6)00 = -ISGeWF^TT 
{'He)ij - 
(^7)00 ^ 



a^(5. 



ij J 



-30G7(7r)ff3^3 ^ 

12G7(7r)ij3^3 ^ q^^^t^hV + 12GT{7T)HHTr^ + \m7{^)H^T:^% 



dt 



a^Sij , 



("^8)00 
i'Hs)ij 



H^ 



, dGs{7:) 
dt 

-,2'^^H^ „ ,,f^l^H^ _ 32^^^i/F 



di 



di2 



(fi 



a ^, 



(A4) 



(A5) 



(A6) 
(A7) 
(A8) 

(A9) 



Appendix B: Explicit forms of £k and {T-ik)ij,y in the static and spherically symmetric space-time 



In this appendix, we give explicit forms of £k and {'Hk)^ii' in the static and spherically symmetric space-time, 
especially in the Schwarzschild space-time. 

In the static and spherically symmetric space-time in ()90p . we find 



^(0) 



'(0) 



(n')^e-^^d^^ (ttO-^c-sa /d$\2 (n')^c-^^ dA d<^ (TrM^e-^A d$ 
7 ^ ^ „ — - -t-5 ^ ^ ^ l^l -41 ^ „ ^^-1-2- ' 



dr'^ 



dr I 



dr dr 



dr 



-20 ^ ^ ^ ^- — —^ — + (third derivative terms) , 



dr 



dr 



10- 



{■n'fe-^^d'^'^ (7r')3e-6A /d$y i-K'fe-^^ dhd^ ^(■K'Y{-K")e-^^d^ 



dr^ 



^r) ''' 



dr dr 



24- 



dr 



U'f^-^^i^ (^/)3-6A^^ „(V)V 

+10^^^^ — -r - 22^^^^ -5- + 2^— ^- 

-|- (third derivative terms) , 

2(-')V^^^ + 2(.')V^^f? '^ 
dr^ V dr 



^, (V)^(V0e-6A JV)3e-6A 



8(.')^e--^^ 4- 12(^^:1!^^ 
ctr dr r dr 



+6(^')(7r")e" ^ - 16 
dr 



^(j) ^^/^2„-4A ^A („l\(„ll\^-iK 



(7r')2e-4^dA , ^ f7r')(7r") 



/N2„-4A 



12- 



.(V)^e 



_/\„-2A 



dr 



/\„-2A 



2(^")e-"^^ - 2(7r')e-"''3^ + 2(7r')e-"'^^ - 4 



dr 



dr 



dr^ 



dr 



{TT')c-^'^dAd^ (7r')e-^A^$ 



dr dr 



dr 



^gK)e2^d$ _ , Jn')e-^^ d-P _ ^ (ttQc-^^ dA _^ 24(^')c"^^ dA 



dr 



dr 



dr 



dr 



(Bl) 



(B2) 

(B3) 
(B4) 
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(B5) 



' 2 2 ' 

y,2 gjj,2 y,2 ^j,2 j,2 \dr J r^ \ dr 

-36^^^ — — + 180^A^ — — + 24 ^ '^ ' —-72 ^ "• ' — . B6 
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Here tt' = — — and tt" = —-ir . Especially for the Schwarzschild metric (j94| , we have 
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In the static and spherically symmetric space-time in (|90p . we obtain 
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